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ABSTRACT 

We analyze observations from the Chandra X-ray Observatory of the south- 
east knot of the Cygnus Loop supernova remnant. In this region, the blast 
wave propagates through an inhomogeneous environment. Extrinsic differences 
and subsequent multiple projections along the line of sight rather than intrinsic 
shock variations, such as fluid instabilities, account for the apparent complex- 
ity of the images. Interactions between the supernova blast wave and density 
enhancements of a large interstellar cloud can produce the morphological and 
spectral characteristics. Most of the X-ray flux arises in such interactions, not in 
the diffuse interior of the supernova remnant. Additional observations at optical 
and radio wavelengths support this account of the existing interstellar medium 
and its role in shaping the Cygnus Loop, and they demonstrate that the south- 
east knot is not a small cloud that the blast wave has engulfed. These data are 
consistent with rapid equilibration of electron and ion temperatures behind the 
shock front, and the current blast wave velocity ~ 330kms~^. Most of this 
area does not show strong evidence for non-equilibrium ionization conditions, 
which may be a consequence of the high densities of the bright emission regions. 

Subject headings: ISM: individual (Cygnus Loop) — shock waves — supernova 
remnants — X-rays: ISM 



1. Introduction 



Supernova remnants play a fundamental role in processing matter and energy in galaxies. 
Supernovae create the heaviest elements, and their remnants mix the enriched ejecta into the 
surrounding interstellar medium (ISM). Supernova remnants (SNRs) are the primary source 
of the hot component of the ISM, and the supernova rate, with the remnants' subsequent 
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evolution, determines whether this component predominates in any particular galaxy. The 
extant environment reciprocally affects the evolution of a SNR. If large, dense clouds populate 
the surroundings with large covering fraction, the blast wave's expansion will stall, limiting 
the amount of thermal energy that is injected into the ISM. 

The Cygnus Loop is a middle-aged (8, 000 years old) supernova remnant, now interact- 
ing with large-scale inhomogeneities in the surrounding ISM. These variations are respon- 
sible for its observable characteristics, including a limb-brightened X-ray shell and strong 
correlations between X-ray and optical surface brightness. Propagating through the low- 
density environment its progenitor evacuated, the undisturbed blast wave has a velocity 
Vs ~ 350 km s^^ and produces low surface brightness X-ray emission with temperature 
kT ^ 0.15 keV {T 1.8 x 10^ K). Upon encountering a large, dense cloud, however, the 
blast wave is decelerated, to Vg < 200kms~^. The post-shock material rapidly cools through 
optical hne emission to T 10^ K. An important consequence of the cloud interaction is the 
development of reflected shocks. These shocks propagate back through the hot, compressed 
interior of the SNR, enhancing X-ray emission (Hester & Cox 1986). 

The northeastern and western limbs of the Cygnus Loop offer two clear examples of 
these interactions between shock fronts and large clouds (Hester, Raymond, & Blair 1994; 
Levenson, Graham, & Walters 2002). Although less prominent than these bright regions, the 
southeast knot is physically similar. Fesen, Kwitter, & Downes (1992) first drew attention to 
this feature, and subsequent work by Graham et al. (1995) demonstrated that it represents 
the early stage of the encounter between the blast wave and a large interstellar cloud. The 
southeast knot is also located near a portion of the primary blast wave traced in Ha at optical 
wavelengths. These are "nonradiative" shocks, which excite Balmer line emission through 
electron collisions in unshocked gas that is predominantly neutral (Chevalier & Raymond 
1978). Thus, they identify sites where the gas is being shocked for the first time. The shocked 
gas has not yet had time to cool radiatively, so the optical spectra lack the emission lines 
that are typical of fully radiative post-shock regions, such as [O III] and [S II], in addition 
to the Balmer lines. 

The Cygnus Loop is nearby and bright, allowing high spatial resolution and high signal- 
to-noise observations to investigate variations and shock evolution on physically-relevant 
scales. At the 440 pc distance of this SNR (Blair et al. 1999), 1" corresponds to a physical 
size of 6.6 x 10^^ cm. Here we analyze new observations of the Cygnus Loop's southeast 
knot obtained with the Chandra X-ray Observatory (Chandra). In addition to superb spa- 
tial resolution, Chandra also provides simultaneous spectroscopic information for direct and 
unambiguous measurement of the distinct physical character of different areas. Combining 
these X-ray data with previous optical and radio measurements, we distinguish temporal 
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evolution of the post-shock medium from external variations in the environment. 

The three-dimensional geometry of the Cygnus Loop comphcates the interpretation of 
observations. Multiple projections of the current blast wave location or cloud interactions can 
emerge along a single line of sight. As a result, areas that appear projected onto the interior 
of the SNR may actually represent newly-shocked material that is located in the foreground 
(or background), not older shocked gas that is genuinely inside the hot SNR interior. In 
the southeast, we illustrate that many of the observed variations are a consequence of these 
projection effects. 

2. Observations and Data Reduction 

We observed the southeast knot of the Cygnus Loop with the Chandra Advanced CCD 
Imaging Spectrometer (ACIS) on 2000 September 1, obtaining spatial resolution around 1" 
over energy E = 0.3-8 keV, and simultaneous spectral resolution E/AE ^ 12. The total 
field of view of the six CCDs is approximately 25' x 17' and covers the brightest X-ray region 
of the southeast knot, some of the diffuse interior, and some slight X-ray enhancements close 
to the blast wave. Figure 1 illustrates these observations in the context of the Cygnus Loop 
as a whole. The bright southeast knot is at the apex of a slight indentation from the near- 
circular boundary of the SNR. Five of the CCDs (10, II, 12, 13, and S2) are front-illuminated 
(FI), having readout electronics that face the incident photons. The sole back-illuminated 
(BI) CCD, S3, is located entirely outside the blast wave boundary. Because the sensitivity 
of the BI and FI CCDs is different, the S3 data are not useful for measuring background 
emission, so we will not discuss them. 

We reprocessed all data from original Level 1 event files using Chandra Interactive Anal- 
ysis of Observations (CIAO) software, version 3.1, removing the 0.5-pixel spatial randomiza- 
tion that is included in standard processing. (See the Chandra Science Center^ for details 
about Chandra data and standard processing procedures.) We applied current calibrations 
(Calibration Database version 2.28), including corrections for charge-transfer inefficiency and 
time-dependent gain variations, to produce the Level 2 event file. We included only good 
events that do not lie on node boundaries, where discrimination of cosmic rays is difficult. 
We examined the lightcurves of background regions and found no significant fiares, so we 
did not reject any additional data from the standard good-time intervals, for a net exposure 
of 47 ks. 
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We created exposure maps to correct soft (0.3-0.6 keV), medium (0.6-0.9 keV), hard 
(0.9-2 keV), and very hard (2-8 kcV) images separately. The total exposure-corrected image 
(Figure 2) includes all these bands, with the raw pixels binned by a factor of two (to 1 
arcsec^) and smoothed by a Gaussian of FWHM=5". Spacecraft dither during the exposure 
is sufficient to cover the small gaps among the four northern CCDs, but a large gap separating 
the southern CCD remains unobserved. 

Spectral variations are apparent in a false-color image constructed from the three softer 
bands (Figure 3). Nearly all the emission above 2 keV is due to background point sources 
unlikely to be related to the SNR, and many of these hard sources arc evident in the blue 
band of the color image. The emission near the blast wave is generally soft, and the low 
surface brightness interior is harder. Stronger hard emission also arises near the site of 
interaction with the dense southeast cloud. 



3. Spectral Modeling 

We select regions for spectral analysis from the color image, identifying extended areas 
of distinct spectral character. Figure 4 shows these locations. The regions are not of uni- 
form size or shape, in order to isolate unique physical characteristics within each aperture. 
The spatially-limited spectral extractions require specific processing and calibration mea- 
surements. As we noted above, we applied corrections for the charge-transfer inefficiency 
and gain variations that have developed over time in creating the Level 2 events file. An- 
other significant change over the course of the Chandra mission is the reduced soft X-ray 
sensitivity, likely the result of build-up of material on the detector. In constructing ancillary 
response files, we used the CIAO contamination correction to account for this variation in 
quantum efficiency. 

Because the Cygnus Loop covers nearly the entire field of view, we measure the back- 
ground emission in comparable observations of the blank sky. The blank sky background 
spectra have been filtered and processed the same way as the data^, and scaled by the rel- 
ative exposure time. A portion of the 13 and S2 detectors covers the exterior of the SNR. 
We compared local background measurements for source spectra within these detectors and 
found no significant difference from the blank-sky backgrounds. For consistency, we use the 
blank-sky background measurements in all cases. 

We fit the spectra in XSPEC (Arnaud 1996), grouping the spectra into bins of a mini- 
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mum of 20 counts, so the errors are normally distributed and statistics are appropriate. 
We use data from 0.5 keV (below which the calibration becomes uncertain) to 2.0 keV, the 
maximum energy of any significant emission. In general, we expect to find thermal emission 
behind the blast wave of this SNR. We use the collisional ionization equilibrium models of 
Borkowski et al. (Hamilton, Chevalier, & Sarazin 1983; Borkowski, Sarazin, & Blondin 1994; 
Liedahl et al. 1995), including dielectronic recombination rates of Mazzotta et al. (1998), 
the VEQUIL model in XSPEC. 

Wc consider abundance variations, which we quote relative to solar values of Anders & 
Grevesse (1989). At the typical temperature of this area of the Cygnus Loop, oxygen and 
iron produce the strongest emission lines. While groups of elements may be expected to 
vary the same way if they are produced through the same modes of stellar evolution — with 
other alpha elements similar to oxygen and other iron group elements similar to iron — other 
members of these nucleosynthetic families are in fact not at all significant in the observed 
spectra. Therefore, we do not include their variation in the spectral modeling. 

Oxygen is the only element whose abundance differs from the solar value, and it is 
consistently underabundant relative to the Anders & Grevesse (1989) value of {0/H)q — 
8.51 X 10^. Our findings are more typical of the local ISM, however. Within 500 pc of 

the Sun, Meyer, Jura, & Cardelli (1998), for example, find gas-phase {0/H) = 3.43 x 10^, 
correcting for updated transition probabilities (Welty et al. 1999, cited in Andre et al. 2003). 

Some spectra demand more complex models to attain adequate fits, requiring additional 
temperature components or non-equilibrium conditions. In the more complex models that 
are presented, the inclusion of additional model components and their free parameters are 
significant at a minimum of the 95% confidence limit, based on an F-test. All quoted errors 
are 90% confidence for one parameter of interest. Figures 5 and 6 show two example spectra 
and their corresponding equilibrium model fits on the same scale. 

3.1. Equilibrium Results 

All single-temperature equilibrium model fits are listed in Table 1, with the values of 
and the number of degrees of freedom (dof) in the last column. This model yields the best 
fits in regions B, D, E, and K. In these regions, the single-temperature equilibrium plasmas 
have temperatures kT = 0.17-0.19 keV and sub-solar oxygen abundance. These regions are 
generally associated with the southeast knot and located near the shocked cloud's optical 
emission. We discuss these spectral results in the context of this inhomogeneous environment 
below (§5), after considering more complex models. 
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Regions G, H, and I require two thermal components, listed in Table 2. The same fore- 
ground column density is applied to both components. These apertures are located in the 
diffuse interior of the SNR, and they are not directly associated with interactions between 
the blast wave and the southeast cloud or the cavity boundary. These areas are noticeably 
bhicr in the energy-coded color image. Each of these regions is fit with a soft component 
similar to the single-temperature regions, as well as a hotter {kT2 ~ 0.6-0.7 keV) compo- 
nent. Qualitatively, we recognize that these apertures cover material that was heated by the 
passage of the original undisturbed blast wave, before it was strongly decelerated in an en- 
counter with higher-density material. At each location, shocked SNR material extends along 
the line of sight, representing a range of original shock velocities. In the two-temperature 
model, however, neither component quantitatively reveals consistent adiabatic evolution of 
the blast wave. We attempted to fit these spectra with a model that explicitly accounts for 
this evolution (the VSEDOV model), but the resulting fits are worse. 



4. Non-Equilibrium Conditions 

4.1. Non-Equilibrium Ionization 

In addition to the equilibrium descriptions above, we also examined non-equilibrium 
conditions that may be relevant in the Cygnus Loop. One important consideration is the 
ionization state relative to the electron temperature. Initially, the post-shock gas is unde- 
rionized compared with its equilibrium value. In the spectrum, the temperature the lines 
indicate, based on ionization state, then differs from the continuum temperature, which the 
electrons determine. The ionization timescale, r = Ugt, parameterizes the degree of equih- 
bration in terms of ric, the electron density, and t, the time elapsed since the shock passage 
(Gorcnstcin, Tucker, & Harnden 1974). The timescale for ionization equilibrium varies for 
each element and is a function of temperature. For example, at T = 2.0 x 10^ K {kT = 0.17 
kcV), the ionization equilibration e-folding timescale for oxygen tieq ~ 3 x lO^^n^^ s (Liedahl 
1999; Shull & van Steenberg 1982). 

Globally, we do not expect the interior of the Cygnus Loop to be in ionization equi- 
librium, with r ~ 3 X 10^°cm~^s for n ~ 0.1 cm~^ and an age of 8,000 years. However, 
the brightest emission is produced in cloud interactions at much higher densities, and are 
therefore more hkely to be equilibrated. At the southeast, the fainter (and therefore lower 
density) edge regions, in which the elapsed time is also shorter, are most likely to be out 
of ionization equilibrium. In the present observations, non-equilibrium planar shock mod- 
els (the VPSHOCK model in XSPEC) improve the spectral fits of regions A, C, F, and J 
compared with the corresponding equilibrium models. In these cases, we include the inte- 
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grated contribution from the time of initial shock (r, = 0), and fit for the current value of 
the ionization timescale. These results are listed in Table 3, and Figure 7 illustrates the 
non-equilibrium ionization example of the region F spectrum. 

Regions A and C are not strongly out of equilibrium, with upper limits of r ~ 5 x 
10^^ cm~^ s. The non-equilibrium effects are stronger in the lower surface brightness (and 
therefore likely lower density) regions F and J. The best-fitting temperature of region J is 
unphysically high, corresponding to an equilibrium shock velocity Vg = 500 km s^^, but it 
is not constrained well. The uncertainties encompass temperatures that are more typical of 
the southeast Cygnus Loop {kT ~ 0.17 keV). 

In an effort to measure changes in the emission behind the shock front and evolution to- 
ward equilibrium, we examined the spectra near the edge of the SNR. The Balmer-dominated 
filaments directly reveal the current (projected) location of the blast wave. The region around 
a = 20^56^40'^, 5 = 30°15'10" offers the clearest signature, with a bright Balmer filament 
marking the shock front (Figure 8). These filaments are not curved, and this region has no 
other indication of comphcating projection effects over nearly 2' toward the SNR's interior, 
where another nonradiative filament is apparent. 

We extracted spectra from four regions located at incrementally increasing distance 
behind the blast wave (L1-L4). The regions are relatively narrow (24" = 1.6 x 10^^ cm) and 
extend 170" parallel to the shock front. We considered fitting the spectra individually and 
together, to explicitly follow the evolution toward equilibrium of the shocked material. In 
fitting the spectra jointly, we constrained the column density, abundance, and equilibrium 
temperature to be the same in all cases. The upper bound on ionization timescale is a free 
parameter for each region, and this value served as the lower bound on r for the adjacent 
interior region. We allowed the normalization to vary independently in the four spectra. 

The non-equilibrium models do not improve the spectral fits in any case. The indepen- 
dent fits tend toward equilibrium, with r > 10^^ cm^'^ s. Furthermore, in both the indepen- 
dent and the combined fits, the ionization parameter does not vary smoothly behind the 
blast wave. These results are surprising because the time since shock passage is minimized 
and the densities are not necessarily extremely high at these edge regions, so non-equilibrium 
conditions would seem likely. The faintness of the easternmost regions may preclude signifi- 
cant measurement of their genuine non-equilibrium conditions. The equilibrium models do 
show systematic temperature variation from the edge to the interior, which we discuss in 
terms of electron-ion temperature disequilibrium or shock deceleration. 
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4.2. Electron-Ion Temperature Disequilibrium 

Immediately following the passage of a strong shock, all particles acquire the same 
velocity distribution. Because their masses are different, electrons and protons are initially 
heated to different temperatures. The temperatures equilibrate through Coulomb collisions 
in the post-shock region. However, plasma instabilities may promote rapid equilibration at 
the shock front (Cargill & Papadopoulos 1988). 

To investigate the role of coUisionless equilibration, we examine the variation of the 
electron temperature behind the shock front in the spectra of regions L1-L4. Although we 
do not measure the ion temperature, it has little direct effect on these spectra. Similar to the 
combined fits above, we simultaneously fit the spectra of regions L1-L4 with the VEQUIL 
model (Figure 9), constraining the column density and abundance to be the same in all 
cases while allowing the temperature and normalization to be free parameters. A single- 
temperature plasma model provides a reasonable fit, with sub-solar oxygen abundance and 
low foreground column density (Table 4). The results are not significantly different from the 
independent fits of these regions (Table 1), but applying the common column density and 
abundance is more appropriate in this interpretation of the evolution of a single shock front. 

Similar to the independent equilibrium fits of these four spectra, we find a gradual 
temperature gradient increasing toward the interior of the SNR. In detail, however, these 
measurements disagree with the predictions of Coulomb equilibration alone. For a simple 
comparison, we use the observed Balmer filament to define the location of the blast wave, 
and the interior kT — 0.17 keV yields the equilibrium electron temperature. In terms of Tj, 
the ion temperature, the electron temperature evolves 

dTe _ Tj - Te 

(Spitzer 1940, cited in Spitzer 1962). The timescale, t^q, is inversely proportional to density, 
and for n = 0.1 cm~"^, t^q ~ 10^^ s as equilibrium Tg = 1.9 x 10® K is approached. Assuming 
the interior region L4 is fully equilibrated, we use its temperature to determine the shock 
velocity (vg ~ 360kms~^), and we then relate the elapsed time to the observable post-shock 
distances. For n — 0.1 cm~^ and Coulomb equilibration alone, we would expect kT — 0.05 
keV near the edge, and a maximum of only 0.09 keV at the innermost location. This 
theoretical prediction sharply contrasts the observed higher temperatures and their more 
gradual rise. 

One reasonable concern is that an equilibrium model is used to determine the electron 
temperature of each of the regions, while nonequilibrium ionization would be expected to 
influence the effective temperatures we measure. Over most of the range of timescales and 
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temperatures we observe, however, the majority of the hne emission is due to O VII transi- 
tions. The absence of higher ionization states does not strongly alter these spectra. Other 
indications of non-equilibrium, such as diagnostic line ratios, cannot be distinguished at the 
present low spectral resolution, so they do not drive the spectral fitting. Furthermore, as 
described above, the non-equilibrium ionization models empirically fail to provide a phys- 
ically reasonable description of the data. Finally and most significantly, the net effect of 
non-equilibrium ionization would be to increase the measured temperatures, requiring yet 
more rapid temperature equilibration. 

Considering the uncertainties in the temperature measurements. Coulomb equilibration 
could be marginally consistent if the density were very high (n = 2cm~^). However, the 
observed flux precludes such a large density. The 0.4 x 3.0 pc^ region LI is offset 26" to the 
interior of the blast wave. The nearly-spherical SNR has a radius of 84', so the line of sight 
through shocked material at the location of LI is 6.8 x 10^^ cm. Prom the observed emission 
of this volume, we find n = 0.1 cm~^. The Coulomb equilibration model assumes constant 
density over the shocked region, yet region LI is most likely the highest density of the four 
areas. It is located at the boundary of the cavity that surrounded the progenitor rather than 
in the evacuated interior. We therefore conclude that Coulomb equilibration is insufficient 
to produce the observed variation of electron temperature as a function of distance behind 
the shock front, and rapid post-shock equilibration is required. 

These results agree with other observations of similar shocks elsewhere in the Cygnus 
Loop, where the initial ion and electron temperatures arc fully equilibrated (Ghavamian et 
al. 2001; Raymond et al. 2003). In general, equilibration appears to decrease with increasing 
shock velocity (Raymond et al. 2003; Rakowski, Ghavamian, & Hughes 2003). The blast wave 
of the Cygnus Loop is relatively slow, in contrast to faster shocks of younger SNRs, which 
do show significant temperature disequilibrium between electrons and ions (e.g., Raymond, 
Blair, & Long 1995; Laming et al. 1996). 

The overall decline in temperature toward the SNR's boundary may be a consequence 
of global blast wave deceleration in the surrounding material. This shell contrasts with the 
rarefied bubble that the progenitor evacuated, through which the blast wave had propagated 
previously. This blast wave deceleration is a genuine dynamical effect, not merely electron- 
ion temperature disequilibrium, and it has been observed around the entire rim of the Cygnus 
Loop (Levenson et al. 1998). Adopting the equilibrium temperature found in region LI at 
the edge yields a current blast wave velocity Vb^u = 330(+12, — 17) kms~^. This value is 
somewhat lower than the shock velocity of the nearby interior region, L4. A single, abrupt 
step of density contrast 1.4 could produce this velocity ratio in the two regions. 
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5. Interpretation 

Data from ROSAT demonstrate that the southeast knot is the resuh of interaction of 
the supernova blast wave with a large-scale interstellar cloud (Graham et al. 1995). Here the 
X-ray rim is indented from the near-circular boundary that can be traced over most of the 
SNR's periphery (Figure 1). Preventing any portion of the projected blast wave edge from 
reaching this circular boundary requires a large obstacle that extends parsecs along the line 
of sight to impede the blast wave. During the early stages of the encounter, X-ray emission 
increases near the interaction site, a conscqiicnce of pressure enhancement when a reflected 
shock propagates back through the previously shocked and compressed SNR interior. If 
the speed of the forward shock propagating through the cloud remains high, the increased 
density in the cloud shock may also enhance X-ray emission. 

The small size of the optical southeast knot mistakenly led Fesen et al. (1992) and 
Klein, McKee, & Colella (1994) to envision a physically small, engulfed cloud. This early 
work also misinterpreted an optical feature (near region B) as the Mach disk that develops 
after the blast wave has passed. As Graham et al. (1995) show, however, the shock front 
propagates through a significantly neutral, and therefore unshocked medium here. This 
analysis — combined with the location of the X-ray enhancement west of the optical emission, 
which requires that the blast wave propagate east — demonstrates that the optical feature 
cannot be the Mach disk. 

These Chandra observations further support the basic interpretation of the southeast 
knot as the early stage of interaction between the blast wave and a large-scale cloud. A small, 
engulfed cloud would not produce the extended and significant X-ray enhancements that are 
observed interior of the optical cloud, and the Mach disk would not be so luminous. Also, 
material stripped from such a cloud would not generate the bright X-rays toward the exterior 
at a similar temperature. We interpret the present results in the context of interactions with 
the large cloud and dense material located at the boundary of the SNR. 

While a large obstacle is present, the X-ray spectral variations indicate that its surface is 
not uniform. The X-ray emission shows significant multiple projections on sub-parsec scales, 
similar to the very small scale structure within the optical knot that high-resolution optical 
images show (Levenson & Graham 2001). Regions B and C are sites at the immediate edges 
the cloud-shock interaction, with increased X-ray emission located adjacent to the optical 
emission of cloud shocks. 

Spectrally, Region C is hotter than nearly all the other regions fit with a single-temperature 
model, independent of whether the equilibrium or non-equilibrium results are adopted. Be- 
cause of the confusion along the line of sight in the diffuse interior (regions H and I, espe- 
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cially), region L4 offers tfie best measurement of the undisturbed blast wave for comparison. 
It is cooler than region C, with kT = 0.17 keV. At region C, the greater intensity com- 
bined with higher temperature is consistent with the hotter contribution due to a reflected 
shock. This characteristic temperature increase would not be observed during the late stage 
of interaction (after the blast wave has completely engulfed a cloud). A second component 
of the early interaction, the forward shock that propagates through the dense cloud, is not 
directly detected in the spectrum, hkely because the shock front is rapidly decelerated below 
X-ray-emitting temperatures. However, this immediate interaction site does include a softer 
component that appears as mottled red patches in the color image. These soft X-rays within 
aperture C contrast with the exclusively hotter material that extends to the west, evident 
as a bright blue/green area. Spectrally we find kT — 0.193 ± 0.005 keV in the western area, 
without much emission below 0.6 keV. 

In region B, we would also expect to find evidence of the hotter reflected shock or, 
if the eastward-moving cloud shock remains fast enough, the continuation of the slower 
forward shock. Surprisingly, the spectrum of region B does not exhibit a high temperature 
similar to region C. It is somewhat but not significantly hotter than region L4. The column 
density measured at B is higher than that of all other regions fit with single-temperature 
equilibrium models. The temperature does increase {kT = 0.181 ± 0.004 keV) if the column 
density is fixed at a value more typical of the surrounding areas (1 x 10^° cm~^), but the fit 
is worse {x^/i^ — 116/69). Alternatively, a weaker hot component {kT — 0.20 keV) may be 
included in addition to the dominant kT = 0.18 keV component and high column density, 
but it does not significantly improve the fit. Without definitive spectral identification, the 
X-ray morphology and the association with optical emission provide the strongest evidence 
that region B represents the interaction between the blast wave and a facet of the large 
interstellar cloud that is obvious at optical wavelengths. The narrow area of high surface 
brightness distinguishes region B from the diffuse SNR interior or a projected edge of the 
undisturbed blast wave. This interaction interpretation contrasts alternative accounts of this 
feature as either the Mach disk of a prior shock or the emission from a dense shell that the 
SNR swept up while propagating through the homogeneous ISM. 

Region D is coincident with the optical cloud, in which shocks are slowed to Vs < 
200kms~^ and therefore do not produce X-rays. This region shows fittle or no additional 
intrinsic absorption, so the X-ray-emitting region is likely located in the foreground of the 
cloud along the same line of sight. The volume of hot gas is therefore smaller than the full 
line of sight through the SNR, resulting in a lower surface brightness. 

The highest surface brightness area we observe is region A, yet a single-temperature 
model fits the spectrum well. (As in region C, the non-equilibrium results are similar and 
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do not strongly rule out equilibrium conditions.) Here multiple projections can account for 
the increased flux, with a greater column of hot gas located along a single line of sight. This 
bright spot is located near the apex of two nonradiative shocks. The Balmer filaments mark 
two separate projections of the blast wave, each appearing when the curved blast wave is 
viewed at a tangency, through a long line of sight, so we expect the X-ray surface brightness 
to double behind them. The observed surface brightness is increased by a factor of two 
relative to the immediate surroundings. Region A is in fact spectrally most similar to region 
K, which offers a clear edge-on view immediately behind a single Balmer filament. The 
surface brightness of region A is approximately 3 times the surface brightness of region K. 

Observations at 21 cm (Leahy 2002) show the distribution of warm atomic hydrogen 
near the southeast knot and also indicate that shocks have passed just beyond the X-ray 
knot. On large (degree) scales, the neutral hydrogen is distributed broadly around the 
outside of the SNR. At the southeast knot, H I is adjacent to the bright X-ray emission, but 
does not overlap it (Figure 10). This morphology suggests that a portion of the projected 
blast wave has just passed the X-ray knot, and the complete three-dimensional extent of the 
blast wave is not located far outside the X-ray knot. Farther north, the atomic hydrogen is 
offset to the east of the SNR blast wave, which the X-rays mark. This separation indicates 
that the neutral surrounding medium existed prior to the supernova, and it is not a shell of 
ambient material the SNR swept up. While the distance to the atomic material is uncertain, 
the radio emission observed in this velocity range over the whole Cygnus Loop is strongly 
correlated with optical and X-ray features that are characteristic of interactions of the blast 
wave with interstellar clouds, which indicates that this observed neutral medium and the 
SNR are genuinely related (Leahy 2002). 

The connection with the surrounding atomic material suggests that the X-ray- and 
optically-bright areas are the tip of the much larger southeast cloud. The cloud, detectable 
in H I, may also have a substantial cold or molecular component. We estimate the lower 
limit on the atomic density very roughly. Over 1° scales, Nhi ~ 2 x lO^'^cm^^. For a line of 
sight depth of 1°, we find uhi ~ 8cm~^. 

The Ha image and Chandra spectroscopy together show the location of many shock 
front projections. We illustrate these schematically in Figure 11. The easternmost Balmer 
filaments are long, arising where the blast wave is relatively undisturbed over large distances. 
The absolute projected edge is significantly indented relative to the near-spherical boundary 
of the Cygnus Loop, so a large obstacle that extends several parsecs along the line of sight 
impedes the blast wave. The filaments in the projected interior of the SNR are shorter and 
produced near smaller-scale inhomogeneities, which do not preserve a tangent view through 
the blast wave over large distances. 
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The Balmer filament near the cloud interaction of region B is offset forward (east) of 
the X-ray enhancement. The X-rays are due in part to the shock reflected at region C, 
which moves west. No Balmcr filament is associated with the cloud shock of region C. The 
X-ray image illustrates that this shock front is extremely curved. Maintaining a large line 
of sight through the edge of the blast wave while it is strongly deflected around an obstacle 
is unlikely. 

6. Conclusions 

These Chandra observations of the southeast knot of the Cygnus Loop reveal details 
of the encounter between the supernova blast wave and a large interstellar cloud. Most of 
the X-ray flux arises in the interaction, with the increased density of the cloud material and 
reflected shocks that further compress and heat previously-shocked gas. The diffuse interior, 
through which the undisturbed blast wave passed, exhibits significant higher-energy emission, 
but having lower density, it provides less of the total fiux within the field of view. 

On this small scale, as throughout this SNR as a whole, the environment determines the 
X-ray appearance of the Cygnus Loop. The complicated morphology is fundamentally a con- 
sequence of geometry, with significant multiple projections along the line of sight. Overall, 
the X-ray spectra are simple and generally indicate coUisional equilibrium in the post-shock 
regions, even when the time elapsed since the passage of the blast wave is short. Only 
two regions show strong evidence for non-equilibrium ionization, and electron-ion tempera- 
ture equilibration behind the shock front is rapid. Optical and radio data further support 
interpretation of the observed spectral variations in terms of extrinsic properties of the en- 
vironment. Thus, in accounting for the complex appearance of the southeastern Cygnus 
Loop, we do not find evidence for complex shock physics, such as fiuid instabilities, or signif- 
icant intrinsic variations in the evolution of the blast wave. Higher resolution spectroscopy 
at similar spatial resolution could yet reveal such complex physics, but these data clearly 
show significant spectral variations on small spatial scales, which we attribute to the SNR's 
environment. 

We thank the referee, P. Slane, for suggestions that significantly improved this work. 
We thank D. Leahy for providing an electronic version of the H I map. This research was 
supported by Chandra grant GOO-1121 and NSF CAREER award AST-0237291. 
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Table 1. Single-Temperature Spectral Models 



Region 










Zo/Zq 


Counts 


xVdof 


A 


1.0 ±0.3 





168 ±0.003 


6 1+^-2 
"•-■--i.o 


'-'•'-''-'-0.04 


5621 


78/74 


B 


1.6 ±0.3 





175 ±0.005 


2 q+O-8 


n 41 +0-08 


2660 


106/68 


n 


n o -1- n Q 
U.y ± U.o 


n 
u 


iyi it U.UUO 


i.o ± U.4 


n p>i +0-08 

"-'•'-'-'--0.07 


O ( oo 


oo /79 

yo/ ( z 


D 


om 





188 ±0.005 


1.7 ±0.4 


'-'•'J"-' -0.05 


3059 


102/69 


E 


1 9+0.5 
-'-•^-0.4 





-■-"^-0.006 


2 Q+0.8 


"-"•^^-0.09 


2029 


56/55 


F 


o.slJi 





-1 4q+o.oo8 

-■-^^-0.015 




"-'•'J' -0.09 


693 


37/22 


G 


2.5 ±0.2 





-1 oo+O.OOS 
-■-o^ -0.004 




'-'■'-'^-0.05 


4981 


173/86 


H 


Q+0.3 
-^■"J-0.2 





189 ± 0.004 


r 7+1.1 

<J-' -0.9 


"-'■^' -0.06 


5410 


188/92 


I 


1.8 ±0.3 





-1 7O+0.003 
'°-0.004 


5 2+1-1 

^•^-0.8 


4S+0-07 

'-'•^'J-0.06 


4935 


188/88 


J 




'J-'J-O.O: 





164+0-004 
^^^-0.007 


1 8+0-9 


n 09+0.08 
^•'-'^-0.05 


2128 


69/51 


K 


1.2 ±0.4 





i«p;+0.004 
-■-"^-0.005 




0.381°:°? 


3095 


62/69 


LI 


1 4+2.6 
-■-•^-0.9: 





12 ±0.03 


1 4+20 
-■-•^-1.4: 


0.38l°:^« 


294 


10/9 


L2 


U-O-O.O: 





-I4C-+0.007 


1 n+140 

-■-•U-o.i 


n o-l +0.09 
"-'•'-'-'--0.08 


653 


27/20 


L3 







1 CI +0.005 
-'-^-'--0.01 


1 4+1-2 

-■-•^-0.1 


0^29l°:°,^ 


1080 


38/34 


L4 


0.8l°i 





ific+o.ooe 

-'-"-''-'-0.008 


-■-■0-0.4 




1724 


72/48 



^Column density, in units of lO^icm"^. 
''Temperature of thermal plasma, in keV. 

^Normalization of thermal component in units of 10~^i^, -where K — 
(10"!"^/ {AnD"^)) f neUndV, D is the distance to the source (cm), Ue is the elec- 
tron density (cm~^), and uh is the hydrogen density (cm~^). 



Note. — Parameters and errors constrained by hard limits are marked -with 
a colon. 



-17- 



Table 2. Two-Temperature Spectral Models 



Region 














Counts 




G 
H 
I 


-■-•^-0.5 


U.i»U_o.o04 
U.i ^ d_0.005 


4.4ili 

2.8;^:^ 


0.48^H^ 
0.531°:°? 

0.401^6^ 


'-'•'-''-'-0.07 

0.69+°:°^ 
0.541°:°? 


8.1^?:^ 
6.911:? 

' -3-1.4 


4981 
5410 
4935 


130/84 
152/90 
120/86 



'^Column density, in units of 10^^ cm~^. 
"^Temperature of thermal plasma, in keV. 

^Normalization of thermal component in units of 1Q~^K, where K — 
(10~^^/(47r£)^)) J neUndV, D is the distance to the source (cm), Ue is the electron 
density (cm"^), and uh is the hydrogen density (cm~^). 



^Normalization of thermal component in units of 10 ^K. 



Table 3. Plane- Parallel Shock Models 



Region 






A'' 


Zq/Zq 




Counts 


xVdof 


A 
C 
F 
J 


1 q+0.5 
l-7lgi 

2.51?:^ 


Q -, 74+0.007 

'^-0.009 

0.198lo:oo9 
0.32l°:°6 


' -^-1.3 
9 1 +0.7 

9 c+180 

0.31^:? 


0.48 ±0.12 
0.67 ±0.2 
0.56 ±0.2 
0.46l°:J7 


1-01^.6 

0.241°:^ 

o.o2ii°:^^5 


5621 
3736 
693 
2128 


74/73 
94/71 
28/21 
63/50 



^Column density, in units of lO^^cm"^. 
''Temperature of thermal plasma, in keV. 

^Normalization of thermal component in units of 10"^-?^, where K = 
(10~^''/(47ri5^)) J neUndV, D is the distance to the source (cm), Ue is the electron den- 
sity (cm~^), and nu is the hydrogen density (cm^'^). 

'^Ionization timescale, in units of 10^^ cm~^ s. 



Note. — 



Parameters and errors constrained by hard limits are marked with a colon. 
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Table 4. Combined Edge Region Spectral Model 



Region Nh^ kT^ Zq/Zq Counts xVdof 

0.35lg:gl 153/117 
294 
653 
1080 
1724 



^Column density, in units of lO^^cm"^. 
''Temperature of thermal plasma, in keV. 

^Normalization of thermal component in units of 1(}~^K, where K = 
(10~^'*/(47rD^)) J neUndV, D is the distance to the source (cm), Ue is 
the electron density (cm~^), and uh is the hydrogen density (cm~^). 

Note. — Parameters and errors constrained by hard limits are marked 
with a colon. 
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LI 0.135^Q 0.5^Q 2 

L2 0.146lHS^ 0.9+°i 

TO n 1 t;Q+0.005 1 q+0.6 

UO U.lOO_Q 008 -'-•<J-0.2 

L4 0.1651HS1 1.6i°| 
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Fig. 1. — ACIS field of view overlaid on mosaic of ROSAT HRI soft X-ray image (Levenson 
et al. 1997, completed with later observations). 
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Fig. 2. — Total (0.3-8 keV) exposure-corrected image, scaled linearly from (white) to 
7 X 10~^ photons cm~^ arcsec"^ (black). The large gap that separates the four northern 
ACIS detectors from the southernmost one is unobserved. 



30^ 00" 

20^ 56"' 30' 00' 55"' 30' 
Right Ascension (J2000) 

Fig. 3. — False-color image, constructed from soft (red; 0.3-0.6 keV), medium (green; 0.6-0.9 
keV), and hard (blue; 0.9-2 keV) exposure-corrected images. The individual images have 
been smoothed by a Gaussian of FWHM=5" and scaled linearly, from to 4, 1, and 0.3 
X 10~^ photons cm^^ arcsec"^, respectively. 
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Fig. 4. — Spectroscopic apertures are identified on the false-color Chandra image. 
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Fig. 5. — ACIS spectrum of region A (crosses) with the equihbrium model (histogram). The 
thermal plasma is at temperature 0.17 keV, absorbed by Nh = 1-0 x 10^^ cm~^. 
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Fig. 6. — ACIS spectrum of region H (crosses) with the best-fitting model (histogram), 
on the same scale as Figure 5. The model includes two thermal equilibrium plasmas, at 
temperatures 0.16 and 0.57 keV, absorbed by Nh = 2.1 x 10^^ cm"^. 
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Fig. 7. — ACIS spectrum of region F (crosses) with the best- fitting model (histogram). The 
non-equihbrium ionization plane-parallel shock model has equilibrium temperature kT = 
0.15 keV, ionization timescale r — 2Ax 10^^ cm~^ s, and is absorbed by Nh — 2.5 x 10^^ cm~^. 
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Right Ascension (J 2000) 

Fig. 8. — Chandra total-band contours overlaid on the Ha image of the southeast knot (from 
Levenson et al. 1998). The optical image is scaled logarithmically, and the X-ray contours 
are linear. Dashed lines mark the boundary of the ACIS field of view. The brightest X-ray 
emission (Region A) is located south of the optically-bright cloud. The X-ray enhancement 
of Region B is coincident with a radiative shock. The narrow Ha filaments toward the 
exterior of the SNR mark the boundary of the blast wave. 
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Fig. 9. — ACIS spectra (crosses) of regions LI through L4. Region LI (blue) is located 
2 X 10^^ cm behind the blast wave. Regions L2, L3, and L4 (green, red, and black, respec- 
tively) are each offset a further 2 x 10^'' cm behind the preceding region. The column density 
and abundance are the same in all regions, while temperature and normalization vary in the 
ionization equilibrium models (histograms). The observed temperature gradient is inconsis- 
tent with Coulomb collisions alone to equilibrate the electron and ion temperatures behind 
the shock front. 
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Fig. 10. — Contours of neutral hydrogen 21 cm emission (from Leahy 2002) overlaid on 
the ROSAT HRI image. Neutral material extends up to the southeast X-ray knot but 
survives immediately outside, indicating that the blast wave has not fully progressed (in 
all projections) far beyond the X-ray knot. Farther north, H I emission is offset east of 
the current blast wave location, which the X-rays indicate. Contour levels are Nhi = 2.0, 
2.3, 2.7, and 3.0 xlO^'^cm^^, in the 2.64-km-s^-'^-widc band centered on v^sR = 3.2 km s^^. 
Only warm, not cold, neutral material produces the 21 cm emission, so these contours arc 
lower limits on the total hydrogen column density. The angular resolution of the 21 cm 
observations is 2' (E-W) x4' (N-S), and they do not cover the area marked with dotted lines. 
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Fig. 11. — Projected shock front locations on the total-band image, scaled linearly from 
(white) to 5 X 10~^ photons cm~^ arcsec"^ (black). Portions of the projected shock front 
that are observed as Balmer filaments are drawn with solid red lines. The dominant high- 
temperature component measured in region C provides evidence for the shock front that 
interacts with the optical cloud; we identify the extended shock front (dashed red line) 
based on the similar total-band and color images. Overlapping projections of multiple shocks 
account for the brightest — yet single-temperature — area, indicating that the observed shock 
fronts continue into the projected interior of the SNR, marked in yellow. 



